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^ ! 

O ■ ABSTRACT 

O ■ 

<N : 

^! We investigate the influence of the inner profile of lens objects on gravita- 

^ ! tional lens statistics taking into account of the effect of magnification bias and 

^ ! both the evolution and the scatter of halo profiles. We take the dark halos as 

CN ! the lens objects and consider the following three models for the density profile of 

^ i dark halos; SIS (singular isothermal sphere), the NFW (Navarro Frenk White) 

^ i profile, and the generalized NFW profile which has a different slope at smaller 

I radii. The mass function of dark halos is assumed to be given by the Press- 

^ I Schechter function. We find that magnification bias for the NFW profile is order 

O ! of magnitude larger than that for SIS. We estimate the sensitivity of the lensing 

^ I probability of distant sources to the inner profile of lenses and to the cosmological 

"Th I parameters. It turns out that the lensing probability is strongly dependent on 
the inner density profile as well as on the cosmological constant. We compare the 

O I predictions with the largest observational sample, the Cosmic Lens All-Sky Sur- 

I vey. The absence or presence of large splitting events in larger surveys currently 

I underway such as the 2dF and SDSS could set constraints on the inner density 
profile of dark halos. 



Subject headings: cosmology: gravitational lensing — dark matter: clusters 



1. Introduction 

It has been known since the 1930s that dark matter is the gravitationally main compo- 
nent in a variety of astrophysical objects such as galaxies and clusters of galaxies, but the 
nature of dark matter still eludes us. ( [Peebles 1993|) . Recently, the systematic discrepancy 



between numerical simulations and observations regarding the inner density profile of dark 
halos has been reported. For the inner density profile of dark halos p(r) oc r""", numerical 
simulations suggest the steeper profile a ~ 1 — 1.5, while observations suggest the shal- 
lower profile a ~ — 1. Numerical simulations of CDM halos by Navarro, Frenk & White 
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(1996,1997, hereafter NFW) have shown that the density profile has the "universal" form 
p oc r~^(r + r^)"^, where is the scale length, irrespective of the cosmological parameters, 
the initial power spectrum and the formation histories. Following NFW, higher-resolution 
simulations have been performed, and the results suggest the steeper inner slope p(r) oc r~^^^ 
( [Moore et al. 19^9i [Fukushige fc Makino 200l|) . On the other hand, observations of rota- 
tion curves of the spiral and the low surface brightness galaxies indicate the shallower halo 
profile a ~ - 1 (|van den Bosch et al 2000| , |de Blok et al. 200 1| , [Borriello fc Salucci 200l| ). 



The mass distribution of the cluster CL0024+1654 is reconstructed from lensed images and 
indicates the flat core ([Tyson et al. 1998|) . It is important to resolve the discrepancy to get 
clues to the nature of dark matter. 

In this regard, strong gravitational lensing effects can provide an important method to 
probe the nature of dark matter. Strong lensing probes the inner dense region of lens objects, 
and the impact parameter is estimated as 

e-6.2/.,o- kpc^^, (1) 

where Di is the distance to the lens and 9 is the image separation. So it is sensitive to the 
inner profile of lenses and can be a useful method for probing the inner profile of dark halos. 
In fact, ^ in Eq.(|l|) is comparable to the scale length in the NFW profile p oc r~^{r + rs)~'^ 
( [Navarro, Frenk fc White 1997[ ). Dark halos is gravitationally dominant in galactic halos 
and clusters of galaxies. However, at the inner region of galaxies, baryonic components 
such as bulge and disk are also gravitationally dominant, and dissipation processes among 
the baryons are important. On the other hand, since clusters of galaxies are formed only 
recently, the baryonic component distributes broadly as gas ( [Rees fc Ostriker 1977[ ) and the 
radial distribution of gas mass is similar to the total mass ( [Einasto fc Einasto 2000|) . Hence, 
in order to study strong gravitational lensing without the need to include the gravitational 
effects of baryonic components, we shall mainly concentrate on clusters of galaxies as the lens 
objects and look for large-separation images (the effect of baryons on lensing was studied in 
( Porciani fc Madau 200q , [Kochanek fc White 200l[ , jKeeton 2001]) ). 



In this paper, we examine the effect of the inner density profile of lens objects on grav- 
itational lens statistics, including the effect of magnification bias and both the evolution 
and the scatter of halo profiles. Statistics of gravitational lensing of QSOs provides a useful 
tool to set constraints on the cosmological constant ([Fukugita et al. 1992| , [Kochanek 1996[ ). 
However, it depends on the lens model such as the profile of lenses and its number density as 
well as cosmology. Hence, in using it as a tool to limit the cosmological constant, we must be 
careful about the uncertainties concerning the lens model. We estimate the sensitivity of the 
lensing probability of distant sources to the inner profile of lenses and to the cosmological 
constant. We consider three kinds of density profile of lens objects : SIS (p(r) oc r~^), the 
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NFW {p{r) oc for the inner profile), the generahzed NFW {p{r) oc r~" for the inner 
profile). Only smooth and spherical models are considered. Subclumps (cluster galaxies) in 
clusters do not affect the cross section of lensing ( [blores et at. 20UI1| ). Nonsphericity can 
affect the relative frequency of four- image lenses ( [Rusin fc 'I'egmark 2001 ). The distribution 



of lenses is taken to be the Press-Schechter function. We compare the predictions with the 
large observational data, CLASS (the Cosmic Lens All-Sky Survey), and predict for larger 
surveys such as 2dF and SDSS. Multiple images of large separation angles are expected to 
be caused by clusters of galaxies, so when we compare the observational data, we will use 
the large angle images {9 > 6"). The statistics of wide-separation lenses (^ > lOkpc) has 
been used to probe the distributions of mass inhomogeneities derived from the CDM sce- 
narios, since one does not need to be concerned with the physics of baryonic components 
and bias ( |Wambsganss et al. 1995| ). By studying the lensing properties in this regime as- 
suming the CDM scenarios, theoretical calculations based on N-body simulations ( |Cen et al. 



1994| , [Wambsganss et al. 1995| , [Wambsganss et al. 1998D and semi-analytical method using 



the Press-Schechter function ( [Narayan fc White 1988|, Kochanek 1995| , Nakamura fc Suto 
|1997| , [Mortlock fc Webster 2000|) have been used to place limits on cosmological parameters. 
Similarly, incorporating realistic input for mass profile and number density of clusters, ana- 



lytical calculations have been used to set constraint to cosmological models (fTomita 1996 ) 



and lens models ( [Maoz et al. 199^ ). We pay particular attention to the influence of the 



inner lens structure on the statistics of wide-separation lenses (for earlier discussion before 
NFW profile, see b'lores fc Primack 1996| ). 



Recently a similar analysis has been performed by several authors QWyithe, Turner"fc 



Spergel 200q , [Fox fc Pen 200T| , |Li fc Ostriker 200q , [Keeton fc Madau 2001|) . Wyithe, Turner 



& Spergel (2001) used the generalized NFW lens model and suggested the optical depth 
to multiple imaging of the distant sources is very sensitive to the inner lens profile, but no 
comparison with observational data was made. Li & Ostriker (2000) found that the lensing 
probability is very sensitive to the density profile of lenses, and somewhat less so to the 
cosmological parameters such as the mean mass density in the universe and the amplitude 
of primordial fluctuations. However, they did not take into account the scatter of dark halo 
profiles (concentration parameter) and the uncertainty in the treatment of magnification bias. 
Keeton & Madau (2001) found that the number of predicted lenses is strongly correlated 
with the core mass fraction. In this paper, we present a systematic study of the effect of 
the inner dark halo profile and cosmological parameters on gravitational lens statistics. It is 
very important to be careful about both the effect of magnification bias, which depends on 
the lens profile and magnification of each images, and both the evolution and the scatter of 
halo profiles in N-body simulations ( [Bullock et al. 2001|) . 

This paper is organized as follows. In section 2, we briefly summarize the basic formulae 
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of gravitational lens statistics for various lens models. In section 3, we examine the effect 
of magnification bias and compare the theoretical prediction with observation. In section 4, 
we estimate the number of lensed images expected in larger surveys. Finally, in section 5, 
we summarize the main results of this paper. We use the units of c = G = 1. 



2. Basics 
2.1. SIS Lens 

The SIS (singular isothermal sphere) model is frequently used in the lensing analysis, 
since it is supported by observed flat rotation curves and moreover the density profile is 
very simple and quantities related to gravitational lensing can be written in simple analytic 
forms ( P?urner, Ostriker fc Gott 1984|) . SIS is characterized by the one-dimensional velocity 
dispersion f , and the density profile is p(r) = f^/27rr^ . The lens equation leads to the two 
solutions (image positions) at x± = y ± 1 with magnifications fj,± = \{y / x±){dy/ dxj-)\~'^ = 
l/y ± 1, where dimensionless quantities x and y are the impact parameter divided by the 
Einstein radius in the lens plane and the source position divided by it in the source plane 
( [Schneider, Ehlers fc I'alco 1992| ). Magnification of the brighter (fainter) image is /i+ (//_), 
and the total magnification is fj,{y) = 2/y. The splitting angle between the two images is 
9 = Snv'^Dji^s/Ds, where D^, Ds and D^s are the angular diameter distances between the 
observer, lens and source. We shall adopt the so-called filled beam distance (pyer fc Roeder 



1973|) , since the ray shooting in an inhomogeneous universe created by N-body simulations 



is consistent with it ([Tomita 1998| , [Tomita, Asada fc Hamana 1999| ) 



When the light ray from the source passes near the lens object, if < 1, double images 
form. We define the cross section cr(f , zl, zs) as the area of a region in the source plane which 
satisfies the following criteria; (i)double images are formed, and (ii)the magnification is larger 
than the minimum amplification /i*. The second condition is needed for the calculation of 
magnification bias (see Eq. ([l5|) below). Then, the cross section is given by, 

a{v, zl, Zs) = 7rr| X 2 / dyyQ{fi{y) - fi^), (2) 

Jo 

where rE{= 4:Ttv^Dls) is the Einstein radius in the source plane, and 0(x) is the step 
function. 
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2.2. Generalized NFW Lens 



In this section we consider the generahzed NFW profile for dark halos ( |Wyithe, 'lurner 
^ Spergel 2UUll|, |Li fc OstrikiFl 



= / .a 73=^' (3) 



where ps and are parameters which depend on the mass of halo M and redshift z, and 
a(0 < a < 2) is a constant. The scale radius is about 10 kpc on a galactic halo scale 
(M ~ W^Mq) and 100 kpc on a cluster scale (M ~ lO^^M©) (|Navarro, I'renk fc White 



|1997|) . Since the impact parameter ^ is comparable to the scale radius (see Eq.(P), we 
could set a strong constraint on inner slope a in Eq. @ by using strong gravitational lensing. 
We consider the case of a = 0.5, 1, 1.5. 

The concentration parameter c and the characteristic density Sc are respectively defined 

by 

c = (4) 
rs 

r Ps ^vir C , , 

" = 3 /o^rfxx2-(l+x)-3- ^> 

c represents degree of the mass concentration at the inner region of a halo. The virial radius 
of the halo r^r is related to the mass M and redshift z as 

/for.., = 5.643 xlO-h'/3n-VS(i^.,)-i('^j _M_^ _ (g. 



0. 



where A^^^ is the overdensity of the halo {Amr = ISvr^ for EdS {Qq = 1) model, and for 
other cosmological models we use the fitting formulae given in Kitayama fc Suto 1996 and 
Nakamura fc Suto 1997| ). 



For the NFW profile (a = 1), c is fitted by the numerical simulation ([Navarro, Frenk 



White 199?! ). For the generalized NFW profile, following the recent work ([Bullock et al. 



2001] , [Keeton fc Madau 200 1|) , we define c to be c = r„jr/r_2, where r_2 is the radius at 



which the logarithmic slope of the density profile is —2 (r_2 is equivalent to for the NFW 
profile). Thus we obtain c = (2 — a.)~^rmrlfs- 

However, there is a scatter in the concentration parameter c at fixed mass and redshift 
( [Bullock et al. 200Ti |Jing 2000| ). Keeton & Madau (2001) pointed out the importance of 



the scatter of c on the lensing statistics. We adopt a log-normal function as the probability 
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distribution function of c, 

p{c)dc 



exp 



(In C - In Cmed) 

2al 



(ilnc, 



(7) 



with = 0.18 ( IBullock et al. 200 1| , |Jing 2000|) . 



Using the toy model of a Cmed — M relation found by Bullock et al. (2001), we estimate 
the concentration of the halos, 



10 / M 



^med 



l + z VM, 



where the fitting parameters (M*, /3) are (7.0 x lO^^M©, 0.16) for EdS {h = 0.5, fio = 1, Aq = 
0,as = 0.67) model, (2.1 x 10^3^^^0,0.14) for A (h = 0.7, fio = 0.3, Aq = 0.7, cxg = 1) model 
and (2.6 x 1O^3M0,O.19) for open {h = 0.7, Qq = 0.3, Aq = 0,^8 = 0.85) model. Recent 
high resolution numerical simulations ( [Bullock et al. 2001 ) show that c does depend on 
z contrary to the earlier suspicious that c does not vary much with the redshift ([Navarro 
[Frenk fc White 1997|) . 



The lens equation for the halo with the generalized NFW profile is 



y = X — A 



X 



(9) 



where x = ^/vg is the impact parameter in the lens plane), y 
source position in the source plane), and 



A 



dz dO — , 

(cos 6* + ^)3-° ' 

DlDls 
167r ps — r,. 



{DjDs){7]/r,) {r] is 

(10) 
(11) 



For a = 0,1,2, the integration of the Eq.(|TD|) can be carried out analytically (the lens 
equation for the NFW lens was obtained by Bartelmann (1996)). In Fig. 2, we show the lens 
equation for the NFW lens with A = 10. If \y\ < ycriu three images Xi {i = 1,2,3,X3 > 
^2 ^ Xi) form. The image positions Xi = Xi{y) with magnifications Pi{y) can be obtained 
numerically. The image separation angle 6 is defined as the separation between the outer 
two images and depends on the source position y. We use the averaged value in the source 
plane. 



e{M,ZL) 



Vcrit 



d{ny^){x3 - xi 



rs{M,ZL) 

Dr 



(12) 
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The total magnification is the sum of tlie magnifications of tliree images (J2i=i l^i) ^^nd tlie 
magnification of tlie fainter image is smaller magnification of outer two images (min { fii, ^3 
}). Following the case of SIS (see Eq.(0)), the cross section is defined by 

/ Dq \ ^ rVcrit 

a(M, ZL, zs) = n (^-^r,j x 2 dyy e(/i(y) - /i,). (13) 



2.3. Lensing Probability 



The lensing probability for a source at redshift zs is ( Schneider, Ehlers fc Falco 1992|) 

Jo Jo iJc 



-2 -TrT^^L a{M,ZL,zs) Nm{M,zl) 



(14) 



where H{z) = Hq [ + z)^ + {1 — Qq — Ao)(l + z)"^ + Aq ]^''^ is the Hubble parameter at 
redshift z, and Nm{M, zl) is the comoving number density of lenses and is assumed to be 
given by the Press-Schechter mass function. The one-dimensional velocity dispersion v is 

1 /2 

related to the mass through v = {M/2rmr) , and this relates the PS mass function to the 
SIS lens profile. The probability of the image separation angle is obtained by computing 
dP/de. 

Since gravitational lensing causes a magnification of images, lensed sources are over- 
represented in a magnitude-limited sample and the actual lensing probability is enhanced. 
This selection effect is called magnification bias. Let $5(2:5, L) dL be the luminosity function 
of sources. The observed flux S for a lensed source is related to the luminosity L, L = 
4:7i{l + zs)^ -D| {1 + zs)'^~^S, where the factor {1 + zs)'^~^ is the K-correction, which assumes 
that the energy spectrum of source is of the form E oc When one searches for lensed 
source of the observed flux S, the lensing probability increases as 



1 



d 



P{zs] 



(15) 



Including magnification bias, the lensing probability P^{zs, L) for the generalized NFW 
lens is expressed as Eq.(|l4D with a{M, zl, zs) replaced by 

a''{M,ZL,zs,L) = 71 (jfrA x — -/ dy y <^s{zs, L/ fiiy))/ M- (16) 

\Dl ) 9s{zs,L)Jo 

fi{y) can be the magnification of the total images or the fainter image. We will discuss each 
case in Sec. 3.1. 
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2.4. Press-Schechter Function 



We use the Press-Schechter (PS) function ( [Press fc Schechter 1974| ) for computing the 
number density of the halos. The PS mass function is given by 



NMiM,z) dM 



TT M (t{R) 



din a 



dlnM 



exp 



2a\R) 



dM 



(17) 



Here po is the mean mass density of the universe at present and (t{R) is the hnear density 
fluctuation presently on the comoving scale R, 

1 /"oo / 9M \^^^ 

o\R) = — dk P{k) W\kR) , 7^ , (18) 

where Piki) and W{kR) are the power spectrum at present ( [Bardeen tt al. 1986| , [Sugiyama 



1995 ) and the top-hat window function. We normahze oiRu) so that a[R = 8/i~^Mpc) = ag, 
and the critical density contrast is 6crit{z) = l.Q86/Di{z), where Di{z) is the hnear growth 
rate normalized to unity at 2; = 0. 



3. Results 

In order to explore the dependence of the lensing probability on cosmological models, 
we consider three representative models; EdS model {h = 0.5, Qq = 1, Aq = 0, as = 0.67), A 
model {h = 0.7, fio = 0.3, Aq = 0.7, ag = 1), open model {h = 0.7,^0 = 0.3, Aq = 0, o"8 = 
0.85). 



3.1. The Effect of Magnification Bias 

We demonstrate the effect of magnification bias by calculating the biased lensing prob- 
ability Eq.(|^, |15D. In calculating magnification bias, we use an optical QSO luminosity 
function. We adopt the 2dF QSO redshift survey data which include about 6000 QSOs with 
the redshift distribution 0.35 < z < 2.3 ( [Boyle et al. 2000| ). The QSO luminosity function 
is fitted by the two-power-law model 

L,{zs) = LoJO'^i^s+'^^^i. (20) 
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The fitting parameters for A model are given by 

{ci,C2,MB„h,k2) = (3.41, 1.58, -21.14 + 51og/i, 1.36, -0.27), (21) 
$^ = 2.88 X 10~*^/i^Mpc-^mag-\ (22) 

wfiere Mb* is tlie absolute B-band magnitude corresponding to Lq* (Eq.(pOD). We take the 
absolute magnitude of a source is Mb = —25.8 mag in Eq.(^Dp and consider the image 
separation range 6 > 0.3". For the energy spectrum index of source (see section 2.3), we 
use 7 = 0.5 for the optical QSOs ([Boyle et at. 1985|). In Fig. 2, the effect of magnification 



bias for SIS and the NFW is shown for A model. The vertical axis is the lensing probability 
with magnification bias divided by that without it. As the source redshift is higher, the 
amplitude of magnification bias is smaller. This is because the number of fainter QSOs (its 
luminosity is L < L*) is smaller at higher redshift (Eq.(|l^)), so the integration of luminosity 
function Eq.(|16D takes a smaller value. From Fig. 2, we find that the magnification bias effect 
for the NFW is order of magnitude larger than that for SIS. This is due to the fact that 
the magnification for the NFW is divergent aX y = ycrit and y = (see Fig.2). Hence, if 
we attempt to predict lensing frequencies by using the NFW lens model, magnification bias 
should not be ignored. 

In Fig.2, we also compare the case when /i in Eq.(|l^) is the magnification of the total 
images with the case of the fainter image. Depending on the properties the gravitational 
lensing configuration, we should use fi in the bias factor as the magnification of the total 
images (we call "/i total") or the fainter image among the outer two images ("/x fainter"). If 
individual sources in a sample are not examined closely enough to determine whether they 
are lensed or not, the magnification of the fainter image should be used ( [Sasaki fc 'l'akahara| 
1993 , Cen et al. 1994|) because the fainter image should be bright enough to be recognized 



as one of the multiple images. On the contrary, if one searches for lensed source of small 
separation angles, then the total magnification may be relevant, because it is likely that the 
brightness of a lensed source with a small separation is recognized as the total brightness of 
all the images. So as far as images of large separation angles are concerned, it may be better 
to use /i as the magnification of the fainter image. 

From Fig.2, different choice of /i in the bias factor greatly changes the amplitude of 
magnification bias for SIS, but does not change it so much for the NFW; the difference is 
only factor of three. This is because the magnification for the NFW is divergent at y = 
more strongly than aX y = ycrit and the magnifications of both of the outer two images 
(including the fainter image) are divergent at ?/ = (see Fig.l). In the following we will 
consider both cases (/i total or /i fainter), since the degree of the magnification bias depends 
on the method of gravitational lensing search and thus the expected number of lensed sources 
will be in between. 
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3.2. Comparison with the CLASS Data 

In this section, we compare predicted lens statistics with a weU-defined observational 
sample. The Cosmic Lens All-Sky Survey (CLASS) is the largest statistically homogeneous 
search for gravitational lenses ( Browne fc Myers 2000 ). The sample comprises 10,499 flat- 



spectrum radio sources with flux 5* > 30 mJy at 5 GHz, and includes 18 gravitational lenses 
with image separations 0.3" < 9 < 6.0". An explicit search for lenses with image separations 
6.0" < e < 15.0" has found no lenses (|Phillips et al. 2000|) . The flux distribution of flat- 



spectrum radio sources in the CLASS samples can be described as a power-law ( [Rusin fc 
[legmark 200^ 



^sizs,L)dLoc L'^-^dL. (23) 

It is a steeper number-flux relation than predicted by the Dunlop & Peacock (1990) lumi- 
nosity function (whose slope is —1.8 for the faint sources). The redshift distribution of the 
full CLASS sample is not known. Marlow et al. (2000) reported the redshifts for a small 
subsample of 27 sources. We assume that the redshift distribution of the full sample is 
identical to that of the subsample. We need not consider the magnitude distribution of the 
sample, since the lensing probability does not depend on the magnitude of the source for the 
power-law luminosity function. 

In Fig. 3-5, we show the predicted image separation distribution of expected number 
of lensed source's in the CLASS with magniflcation bias for each case of fi. The angular 
resolution of parent survey in the CLASS is very low. The survey flux encompasses all 
the flux of even the widest separation lenses observed to date. The CLASS survey then 
reimages the systems looking for multiple imaging ( [Myers et al. 1995| , [IVlortlock fc Webster 



2000| , |Helbig 2000| ). Hence, for the lensed source, /i is given by summing the fluxes of all 
images (i.e. "/z total" is appropriate). In these flgures, we also show the case of "/i fainter " 
for comparison. We note that the lenses of large separation angle {6 > 6") is expected to be 
lensed by clusters of galaxies (or equally dark halos), so we should compare the theoretical 
prediction with large image separation side in CLASS data. In Fig.3, we take the NEW 
proflle and compare the dependence on the cosmological models. The lensing probability is 
the highest for EdS model. This is due to the fact that the PS mass function is proportional 
to Qq (see Eq.(|l7|)) and the concentration parameter c is the highest (see Eq.(|)). In Eig.4, 
we show the lens model dependence for A model. Since there are no lenses for large image 
separations (6" < 9 < 15") in CLASS sample, the steeper inner proflle (a > 1.5) seems 
disfavored. However, it is preliminary, since we do not know whether the subsample by 
Marlow et al. is a fair sample. 

We estimate the sensitivity of the lensing probability to the model parameters. Using 
Eq. (|T^|T5|) and the CLASS data with the magniflcation of total images, it is estimated around 
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a = 1, (jg = 1, Ao = 1 — f^o = 0.7, w = —1 



5No 



-7.4— + 5.7^^+4.3— -5.7-^ + 0.35 — , (24) 
-7.3 — + 6.3 ^^ + 6.1 —-6.3 —^ + 0.27 — , (25) 

g=l2" ^ Cmed (^8 Aq W 



where w is the equation of state of dark energy {w = —1 for the cosmological constant) and a 
flat FRW model is assumed. Eqs.(p4D and (pS]) indicate clearly that the lensing probability is 
very sensitive to the lens model parameters {a, Cmed) as well as the cosmological parameters 
(Ao,o"8), but not sensitive to dark energy parameter (w). For example, in order to put 
constraint on Aq within (9(10)% accuracy, one needs to determine the inner profile a and the 
concentration parameter c with similar accuracy. The dispersion of concentration parameter 
c is about 0.2 ( |Jing 2000t Pullock et al. 2001|) . However, the current uncertainty in a 
is (9(50)% (a ~ 1 — 1.5). Hence, it may be more useful to use the number of large image 
separation lenses to constrain the inner density profile a. Cosmological parameters Aq and cxg 
could be determined within C(10)% by using CMB, SNIa and number count of clusters data 
( |de Bernardis et al. 2001| , [t'an fc Chiueh 200 1| ). The sensitivity of the lensing probability 
to the model parameters was also estimated by Li & Ostriker (2000), but they assumed a 
single source at 2; = 1.5 and did not include the effect of magnification bias. So the detailed 
comparison may not be so meaningful. However, we note that the dependence on c and the 
cosmological constant parameter is slightly larger than that found by Li & Ostriker (2000). 

In Fig. 5, we show the effect of the scatter in N-body simulation on the image separation 
distribution. The dispersion of Ng, aj, is calculated by using the probability distribution 
function of c (Eq.(|^). We use the NFW profile for A model and compare the amplitude of 
the square of the dispersion ag with the predicted number of lenses Ng. We find that the 
amplitude of ag is comparable to or larger than Ng, so the scatter of halo profiles strongly 
affect the lensing probability. When compared with Fig. 3, we also find that the scatter is 
too large to distinguish the different cosmological models, even if the lens model (a) is fixed. 



4. Prediction for Future Survey 

With the current data, the constraints on the parameters are not sufficient. However, 
we expect that larger surveys currently underway such as the 2dF and SDSS detect a larger 
number of lenses. For example, the Sloan Digital Sky Survey (SDSS) plans a spectroscopic 
survey of 10^ QSOs over vr steradian brighter than i' ~ 19 at 2 < 3.0; at redshift between 
3.0 and about 5.2, the limiting magnitude will be i' ~ 20 ( [York et al. 2000]) . In this section, 
we will make predictions for the SDSS. 
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Let NgdO be the expected number of lensed QSOs with image separation 6 ^ 6 + dO 
within sohd angle tt in the sky. We use the QSOs luminosity function $5 for ^ > 3 in SDSS 
data ( [ban et at. 2001| ), since $5 from 2dF redshift survey is known only for lower redshift 
QSOs {z < 3) ( Iboyle et al. 2UUU| ). $s is fitted by a power-law, 



^Sizs, M1450) = ^llO-0-HMi,,o+26-a{zs-mP+l)}^ ^26) 

where M1450 is the absolute AB magnitude of the quasar continuum at 1450A in the rest 
frame. We assume rrii' = mu^o + 0.7. The fitting parameters for A model are given by ( [I'anI 
et al. 20011) 

($* , a, p) = (2.6 X 10~^/i3Mpc-Wg-\ 0.75, -2.58). (27) 



Similarly, for lower redshift sources [z < 3) the QSO luminosity function (Eq. (p^ , |20|) ) is used 



( Poyle et al. 20001) . Then, using the QSO luminosity function $5, Ng{9) can be calculated 

as 

Noie) = / dz^ dL^ie,z,L)<!>siz,L), (28) 

JO dz JLumiz) dU 

dK ^ ^ il + zfDjzy 

dz ^ H{z) 

where Zmax = 5.2 and Lum^z) is calculated from the limiting magnitude. In Fig.6, we show 
the predicted image separation distribution for SDSS. We use various lens model for A model. 
In this model, the number of QSOs is expected to be about 26,000. In Table 1, the expected 
number of large image separation lenses (6" < 9 < 30") is shown. We find that the ambiguity 
resulting from the treatment of magnification bias is not so large. From this table, we expect 
that the future SDSS data could set constraint on the inner density profile. 



5. Summary and Discussion 

We have examined the influence of the inner density profile of lenses on gravitational lens 
statistics carefully taking into account of the effect of magnification bias and the evolution 
and the scatter in halo profiles. We have estimated the sensitivity of the lensing probability 
to the inner density profile of lenses and to the cosmological constant. We have found 
that lensing probability is strongly dependent on the inner density profile as well as on the 
cosmological constant. We have also shown that magnification bias for the NFW is order of 
magnitude larger than that for SIS. There is an uncertainty in the treatment of magnification 
bias: fainter image should be detected in the survey, or the light from both the fainter and 
brighter images is initially unresolved in a single image and thus the total image should be 
detected. However, for NFW profile, difference between the magnification bias of the fainter 
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image and that of the total image is found to be only by factor of three. In any case, we 
should be careful about magnification bias which strongly depends on the lens profile. We 
have compared the predictions with the CLASS data and suggested that the steeper inner 
profile {a > 1.5) seems disfavored. The absence or presence of large splitting events in larger 
surveys currently underway such as the 2dF and SDSS could set constraints on the inner 
density profile of dark halos. 

Recently, using the arc statistics of gravitational lensing, various authors have examined 
the inner profile of dark halos ( [Bartelmann et al. 1998| , [Molikawa fc Hattori 2000| , [Oguri 



[laruya fc Suto 200 1| ). Comparing with the existing observational data, Molikawa & Hattori 
(2001) and Oguri, Taruya & Suto (2001) suggested that the steeper inner profile of dark 
halos (a > 1 or even a > 1.5) is favored. On the other hand, the absence of large images 
separations in QSOs multiple images in the CLASS sample constrains the inner profile and 
rather disfavors the steeper profile. Combining the arc statistics and the statistics of QSOs 
multiple image, we could narrow the allowed range of the inner profile of dark halos, or more 
interestingly both methods might exhibit discrepancy. 

In any event, larger surveys will produce a lot of QSOs multiple images in near future 
and theoretical development especially concerning the uncertainties of various models will 
be expected. So we will get clues to the nature of dark matter. 

We would like to thank Professor Yasushi Suto, Dr. Ryoichi Nishi for useful discussion 
and comments. One of the authors (RT) also thanks Dr. Atsushi Taruya and Dr. Masamune 
Oguri for useful comments. This work was supported in part by a Grant-in-Aid for Scientific 
Research (No. 13740154) from the Japan Society for the Promotion of Science. 
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SIS 


a = 1.5 


a — 1 


a - 0.5 


Nei6" <9< 30") /x total 


197.5 


31.0 


7.3 


1.6 


/I fainter 


37.5 


10.9 


2.4 


0.43 



Tabic 1: The expected number of large image separation lenses (6" < < 30") for SDSS 
data for A model with magnification bias of the fainter image (bottom) and the total images 
(top). The total QSOs number is expected to be about 26,000. 
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Fig. 1. — The lens equation for the NFW. The case with A = 10 is shown. The horizontal 
axis is X, which is the impact parameter normalized by a scale radius in the lens plane; the 
vertical axis is y, which is the source position normalized by a scale radius in the source 
plane. Multiple images are formed when |y| < Ucrit- 
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Fig. 2. — The amplitude of magnification bias for A model {h = 0.7, Qq = 0.3, Aq = 0.7, as = 
1.0) with the absolute magnitude of a source being Mb = —25.8 mag. The horizontal axis 
is Zs, which is the source redshift; the vertical axis is a ratio lensing probability with the 
magnification bias to that without it. Image separation range 6 > 0.3". The solid (dotted) 
line is the NFW profile for the case of the magnification of the total images (the fainter 
image). The short (long dashed) hne is SIS for the case of the magnification of the total 
images (the fainter image). 



-20- 




Fig. 3. — The distribution of image separations for the NFW lens modeL The assumed 
cosmologies are EdS model(short-dashed: h — 0.5, = l^-^o = 0, (jg = 0.67), A model 
(solid: h — 0.7, — 0.3, Aq = 0.7, as — 1) and open model(long-dashed: h — 0.7, Qq — 
0.3, Ao = 0, (Tg = 0.85). The left figure is for the selection condition in which fainter image 
should be detected in the survey, and the right figure is for the selection condition that the 
light from both the fainter and brighter images is initially unresolved in a single image, "/x 
total" is appropriate in the CLASS. The observational data from CLASS are shown by the 
histogram. 
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Fig. 4. — Same as Fig.3, but for various lens models (dotted: SIS, long-dashed: the gener- 
alized NEW with a — 1.5, solid: NFW {a — 1), short-dashed: the generalized NFW with 
a — 0.5). A model is assumed for cosmology. The observational data from CLASS are shown 
by the histogram. 
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Fig. 5. — The dispersion of the predicted number of lenses caused by the scatter of halo 
profiles in N-body simulation. The solid line is the square of the dispersion (standard devi- 
ation), the dashed line is the averaged distribution of image separations. A model and the 
NFW profile {a — 1) are assumed. 
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Fig. 6. — The image separation distribution expected for SDSS data for various lens models. 
A model is assumed for cosmology. The total QSOs number is estimated to be 26,000. 



